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ABSTRACT 

The second Red-sequence Cluster Survey (RCS-2) is a ~1000 square degree, multi-color imaging 
survey using the square-degree imager, MegaCam, on the Canada- France-Hawaii Telescope (CFHT). 
It is designed to detect clusters of galaxies over the redshift range 0.1 <^ z <^\. The primary aim is 
to build a statistically complete, large 10'*) sample of clusters, covering a sufficiently long redshift 
baseline to be able to place constraints on cosmological parameters via the evolution of the cluster 
mass function. Other main science goals include building a large sample of high surface brightness, 
strongly gravitationally-lensed arcs associated with these clusters, and an unprecedented sample of 
several tens of thousands of galaxy clusters and groups, spanning a large range of halo mass, with 
which to study the properties and evolution of their member galaxies. 

This paper describes the design of the survey and the methodology for acquiring, reducing and 
calibrating the data for the production of high-precision photometric catalogs. We describe the method 
for calibrating our griz imaging data using the colors of the stellar locus and overlapping Two-Micron 
All-Sky Survey (2MASS) photometry. This yields an absolute accuracy of < 0.03 mag on any color and 
«0.05 mag in the r-band magnitude, verified with respect to the Sloan Digital Sky Survey (SDSS). 
Our astrometric calibration is accurate to ^ 0.3" from comparison with SDSS positions. RCS-2 
reaches average Scr point source limiting magnitudes oi griz = [24.4,24.3,23.7,22.8], approximately 
1-2 magnitudes deeper than the SDSS. Due to the queue-scheduled nature of the observations, the 
data are highly uniform and taken in excellent seeing, mostly FWHM^O.7" in the r-band. In addition 
to the main science goals just described, these data form the basis for a number of other planned 
and ongoing projects (including the WiggleZ survey), making RCS-2 an important next-generation 
imaging survey. 

Subject headings: surveys — techniques: photometric — galaxies: clusters: general — galaxies: general 
— cosmology: observations 



1. INTRODUCTION 

Clusters of galaxies are ideal tracers of the largest den- 
sity fluctuations in the Universe, and their abundance 
(and its evolution with cosmic time) may be used to 
place constraints on cosmology (e.g.. Eke et al. 1996). 
They also provide ideal laboratories for studying galaxy 
evolution. Originally used in this way since they con- 
tain a large number of galaxies all in the same location, 
it has since become clear that the properties of their 
member galaxies are markedly different from galaxies in 
the general field (e.g.. Dressier 1980, Balogh et al. 1999, 
Ellingson et al. 2001), implying that mechanisms which 
truncate star formation and transform galaxy morphol- 
ogy operate on cluster scales (e.g., Treu et al. 2003 and 
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references therein). 

Constructing large, well-defined samples of galaxy clus- 
ters has a long and varied history. The first system- 
atic searches involved visual identification of overden- 
sities of optical galaxies on photographic plates (Abell 
1958; Abell et al. 1989). In the 1970s, with the advent 
of X-ray telescopes above the Earth's atmosphere, selec- 
tion of clusters from their extended X-ray emission found 
favour (Mitchell et al. 1976; Serlemitsos et al. 1977). Re- 
cently, a combination of large format CCD detectors and 
objective algorithms to search efficiently for signatures of 
galaxy clusters has led to a revival in the use of optical 
selection in cluster surveys (Postman et al. 1996; Kepner 
et al. 1999; Gal et al. 2000; Gladders & Yee 2000; Gilbank 
et al. 2004) . A variety of techniques have been suggested 
to exploit the expected luminosity and/or color distribu- 
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tion of galaxies in clusters. The big advantage of these 
surveys compared with the older visual searches is that 
the detection method could be automated and charac- 
terised, meaning that the survey selection function could 
be quantified. Arguably the most efficient method is that 
of Gladders & Yee (2000) which uses the fact that the 
cores of galaxy clusters are dominated by galaxies with 
old stellar populations, forming a tight red-sequence in 
color magnitude space (Visvanathan 1978; Bower et al. 
1992). A number of other realisations of red-sequence 
based cluster finding algorithms exist (e.g., Koester et al. 
2007) differing in some details but all relying on accu- 
rate colors from imaging in two or more filters. The 
observed color of this sequence provides an accurate dis- 
tance estimate. The application of this method led to 
the construction of the first Red-sequence Cluster Sur- 
vey (RCS-1, Gladders & Yee 2005), a 72 square degree 
imaging survey in two bands {Rc and z') designed to 
locate galaxy clusters from 0.2 ^ 2; ^ 1.1 using the tech- 
nique of Gladders & Yee (2000). 

Not only is optical selection of galaxy clusters under- 
going a revival, but astronomy in general is entering an 
era of 'survey science' where an unprecedented number of 
wide-field optical (and NIR) survcivs arc currently under- 
way or planned, such as LSST (LSST Science Collabo- 
rations et al. 2009); Pan-STARRS^; UKIDSS (Lawrence 
et al. 2007); and DES^. In addition, many of these wide- 
field optical surveys are specifically targeted at areas 
surveyed for clusters using other methods, such as the 
Blanco Cosmology Survey (High et al. 2010) of the South 
Pole Telescope (SPT, Carlstrom et al. 2009) Sunyaev 
Zel'dovich effect (SZ)-selected cluster survey. For surveys 
using other methods (such as SZ selection), the optical 
data arc critical for the verification of the cluster can- 
didates found and for the determination of photometric 
redshifts. Furthermore, surveying the same areas with 
multiple techniques allows important comparisons of the 
different selection effects and the resulting properties of 
the clusters found (e.g., Donahue et al. 2001; Gilbank 
et al. 2004; Rasmussen et al. 2006). 

In this paper we describe the second Red-sequence 
Cluster Survey (RCS-2), the largest survey of this new 
generation for which imaging has already been com- 
pleted. This builds on the methodology of RCS-1. The 
RCS collaboration has invested a large amount of work 
in attempting to characterize the selection function and 
the properties of clusters selected with this technique. 
Many of these results are directly applicable to RCS-2 
(such as mass-richness calibrations) and so it is useful to 
summarize some of the RCS work to date. 

The efficiency of the selection method employed by the 
red-sequence surveys is that it can locate and estimate 
the redshifts of clusters using only one color (two filter) 
data, given the appropriate choice of filters. It is imprac- 
tical to obtain mass estimates from follow-up observa- 
tions of the ~30 000 clusters which will be found in RCS- 
2, so the survey data themselves must be used to produce 
a proxy for cluster mass. Significant, representative sam- 
ples of clusters from RCS-1 have been followed up using 
a variety of mass estimators such as dynamical (Gilbank 

^ see http://pan-starrs.ifa.hawaii.edu/public/ 
^ see https://www.dajkenergysurvey.org/ 



et al. 2007, 2010), X-ray (Hicks et al. 2008), strong 
and weak-lensing (from an ACS snapshot programme, 
PI:Loh; ACS SNe Cosmology project PLPerlmutter) and 
SZ observations. In this way, the relationship between 
our mass proxy (optical richness from the survey data) 
and cluster mass can be understood. 

One of the primary goals of RCS-1 was to place con- 
straints on cosmological parameters (i^M, Gladders 
et al. 2007) via the growth of the cluster mass function. 
This demonstrated for the first time the feasibility of 
such an approach using an optically-selected cluster sam- 
ple. This approach used the measured relation between 
mass and richness, but also showed that meaningful con- 
straints could be obtained using a self-calibration tech- 
nique (Majumdar & Mohr 2004) to estimate the form 
of this relation from the survey data themselves. These 
authors demonstrate that the best constraints are ob- 
tained when accurate mass estimates are available for a 
subsample of clusters within the survey. It is worth em- 
phasising that even if there is significant scatter in the 
relation between mass and the proxy (as we have found 
for optical richness), it is only important that the size 
of the scatter be well understood. With an order-of- 
magnitudc larger survey than RCS-1, it becomes feasible 
to also constrain the equation of state of dark energy, w, 
(Majumdar & Mohr 2004) and this is in part the mo- 
tivation for RCS-2. RCS-1 also produced a significant 
sample of strongly gravitationally lensed arcs around the 
clusters found. The number and redshift distribution of 
these lensing clusters were used to argue about the phys- 
ical properties of the clusters responsible for their lensing 
cross-section and the relevance of such systems to con- 
straining cosmology (Gladders et al. 2003). The identifi- 
cation of such high surface brightness, strongly-lensed 
galaxies is another primary science driver for RCS-2. 
The massive clusters can be used as gravitational tele- 
scopes for studying high- redshift galaxies (e.g., Pettini 
et al. 2000; Wuyts et al. 2010) which would otherwise be 
too faint to observe in detail. The giant arcs can also 
be used as probes of the properties of the cluster lenses 
themselves. 

With a statistical sample of galaxy clusters, such as in 
RCS-1, it is possible to study the properties of their mem- 
ber galaxies (e.g., their luminosity functions and blue 
fractions) by stacking subsamples built from the survey 
data themselves (Gilbank et al. 2008; Loh et al. 2008). 
With the order-of-magnitude larger RCS-2, it becomes 
feasible to measure much weaker trends and push mea- 
surements of cluster galaxies down to much lower over- 
densities. The addition of photometric redshifts (e.g., 
Hsieh et al. 2005) will allow these techniques to be ex- 
tended to the field environment. Such galaxy evolution 
studies will be explored in future work with RCS-2. 

The outline for this paper is as follows. In §2 we give 
an overview of the survey design and observations; §§3 
& 4 deal with CCD pre-processing, reduction, object de- 
tection and photometry; §5 describes the photometric 
calibration via accurate fits to the star colors in our sur- 
vey fields; §6 outlines the procedure for and accuracy 
of the astrometric calibration. §§7 & 8 describe the in- 
corporation of additional data into our primary RCS-2 
catalogs: i-band data which covers a large (~70%) sub- 
sample of the primary g, r, z survey area; and public 
imaging data from the CFHTLS-Wide survey, respec- 
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tively. §9 describes tlie final cleaning of the photometric 
catalogs: stitching into contiguous patches, removing du- 
plicate data between overlapping pointings, and masking 
of artefacts. §10 summarises and describes ongoing and 
future work for the survey. 

2. SURVEY OVERVIEW 

The wide-field imaging capability of the square degree 
imager, MegaCam (Boulade et al. 2003), on the 3.6-m 
CFHT makes it feasible to carry out a survey covering 
a significant fraction of the sky in a modest amount of 
observing time. Coupled with the exquisite seeing condi- 
tions attainable from the summit of Mauna Kea, such a 
survey can achieve impressive depths/resolution within 
this time. 

The main science drivers for RCS-2 dictate a survey 
area of ^1000 deg^ with the ability to detect galaxy 
clusters out to z^l using the method of Gladders & Yee 
(2000). Following the strategy of RCS-1 (Gladders & Yee 
2005), which used Rc- and z-band imaging, r'- and z'- 
band^ filters are chosen to cover the bulk of this redshift 
range. In order to better distinguish lower-redshift clus- 
ters, (/-band is added to accurately measure the colors of 
galaxies at z^O.4, once the 4000A break begins to move 
blueward of the r-filter. The g-band filter also helps to 
identify strongly gravitationally-lensed arcs around clus- 
ters, since the former tend to be relatively blue. Thus, 
the survey comprises three-color, g, r, z, imaging over the 
whole area, i-band imaging is obtained for the majority 
of this area via a data-exchange with the Canada- France 
High-z Quasar Survey (Willott et al. 2005). 

2.1. Observations 

MegaCam comprises 36 CCDs arranged in a four row x 
nine column grid at the prime focus of CFHT. Each chip 
is 2048x4612 pixels with a pixel scale of 0.187", allow- 
ing it to adequately sample the exquisite seeing possible 
from Mauna Kea. The spacing between CCDs is approx- 
imately 13", with larger gaps (?s80") between the up- 
permost and lowermost rows and the other CCDs. This 
gives sky coverage of 0.96x0.94 deg^ (see Fig. 1). 

Observations were carried out in queue-scheduled 
mode on CFHT on runs between semesters 2003A and 
2007B inclusive. PI imaging time was granted through 
requests to Canadian and Taiwanese agencies. 

Exposure times were set to 4, 8, and 6 minutes in g, 
r, and z respectively. In 0.65" seeing, according to the 
MegaCam exposure time calculator, these should corre- 
spond to 5-cr point source limits oi g w25.3, r w24.8, and 
z «22.5. The depth was set in the two reddest bands 
(r and z) by the requirement to reach « M* + 1 red- 
sequence cluster galaxies at z~l. Part way through our 
survey observations, it was found that 0.65" seeing was 
not available on as many nights as required by our pro- 
gram per semester'*. So, a two-tier strategy for r and 
z-band image quality was adopted in which only half the 
survey would be conducted in the 0.55-0.75" bracket and 
half would be conducted in 0.75-0.90". This means that 

^ Hereafter we use the shorthand of omitting the prime notation 
from the MegaCam filter names. 

^ This was somewhat mitigated by the significant improvement 
in Megacam image quality resulting from the L3 lens being replaced 
in the upside-down position in late 2004. 
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Fig. 1. — Example of a catalog produced for one MegaCam point- 
ing, illustrating the layout of the CCDs. Points show r < 24 objects 
classified as galaxies in one example pointing (0133A0). The 9x4 
grid of the individual detectors is clearly visible, as arc the large 
gaps between the top and bottom rows and the others, described in 
the text. Circular regions absent of galaxies show where our masks 
for bright star haloes (described in §9.1) have been applied. 

the lowest richness clusters will not be detected in the 
worse seeing imaging (due to the reduced depth) all the 
way out to z~l, but richer clusters will. Since the former 
arc exponentially more numerous, this should have neg- 
ligible impact on cosmological constraints, g-band imag- 
ing was always performed in the better-seeing bracket 
in order to preserve the low surface brightness require- 
ment for the detection of strongly-lensed arcs. Our final 
measured depths are detailed in §5.3. 

Single exposures were taken at each pointing position, 
i.e., no dithering was used. This is because sufficient 
depth is achieved in only a single exposure of 4-8 min- 
utes (depending on filter), so the overhead associated 
with reading out the CCD (»2 minutes) becomes a sig- 
nificant fraction of the integration time if this is divided 
into two or more exposures. The priority is to cover as 
large an area of sky as possible to these depths in a given 
amount of observing time. With this strategy the sur- 
vey will contain chip gaps, but these can be dealt with 
(for cluster finding, etc.) via geometric correction fac- 
tors (discussed in §9.2). Gaps in the data also appear 
due to the occasional failure of a CCD within the mosaic 
camera.'"' 

Individual pointings are arranged in a grid pattern in 
discrete patches/fields (discussed in the next section). 
Pointings are named with the patch name, followed by a 
letter designating the column, beginning at 'A' running 
east to west. Next a digit specifies a pointing's location 
in the declination direction, beginning at '0', running 
from south to north. Thus, pointings within patch '0133' 
are labelled '0133A0' to '0133F5', from the south-east 
to north-west corners across the patch. Each pointing 

^ In a small number of cases, problems with the camera electron- 
ics have caused half the mosaic to fail and so the queue observers 
have observed the northern and southern half of a single pointing 
in two separate observations with MegaCam operating in this half 
functional mode. In these cases we have processed the 'half point- 
ings' as if they were separate pointings with an 18 chip camera and 
combined them into the final catalog at the stitching stage (§9.2). 
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overlaps with its neighbors by wl arcniin, typicahy. See 
Fig. 2. 

The original observing strategy was to observe a given 
pointing sequentially in all three filters (5, r, z), but 
the overheads associated with filter changing made this 
relatively inefficient and so, after the first semester, the 
queue observers switched to observing several neighbor- 
ing pointings in a single filter and then repeating these 
in the next filter, i-band imaging is typically performed 
many months later for a given pointing. This is poten- 
tially a consideration for projects wishing to use RCS-2 
for time variability studies. 

2.2. Field descriptions 

When conducting a large survey, it is obviously advan- 
tageous to place fields such that they overlap with other 
large surveys at other wavelengths, whenever possible. 
The RCS-2 survey area is divided into discrete patches 

(where 'patch' refers to a contiguous set of MegaCam 
pointings). The patches are divided approximately into 
two categories: those designed to overlap with other sur- 
veys (such as SWIRE. Lonsdale et al. 2003) are typically 
a 6 X 6 grid (i.e., 6 degrees on a side); and those designed 
to be 'blank' pointings typically 9x9. The larger patches 
are designed to be sufficiently large that the clustering of 
galaxy chistcirs may be used to provide additional con- 
straints on cosmology (e.g., Majumdar & Mohr 2004). 
These larger 'blank' pointings are placed in the footprint 
of the UKIDSS (Lawrence et al. 2007) wherever possible, 
so that we will obtain NIR data for our optical imaging, 
useful for deriving stellar masses of galaxies, for exam- 
ple. The patch size is reduced to ~ 6 x 6 where we target 
other surveys, if the field size of the survey with which 
we overlap is smaller than 6x6 deg^. We choose discrete 
patches distributed across the sky in order to improve 
observing efficiency (by distributing fields in RA), and 
to minimise the effects of cosmic variance (although this 
would only be an issue for the rarest systems, such as the 
most massive clusters, since our individual patches are so 
large - the largest patches in RCS-2 are already bigger 
than the entire RCS-1!). Most of our fields are equatorial 
in order to maximize follow-up from telescopes in both 
hemispheres. 

All patches are chosen to be in regions of low Galac- 
tic extinction (see Fig. 3). Another consideration is the 
avoidance of bright stars. Using the Bright Star Cat- 
alogue (HofHeit 1964), patches are shifted slightly to 
minimise the number of stars they contain in the range 
2 < m,v < 4 wherever possible. Stars this bright (and 
indeed somewhat fainter) can cause large reflection arte- 
facts (see §9.1) which can render significant areas of the 
image unusable. Stars brighter than mv < 2 are po- 
tentially hazardous to the detectors and so pointings in 
our grid which would overlap with such stars are not ob- 
served. 

A list of the patches observed is given in Table 1. Patch 
1040 covers the Lockman Hole, also surveyed by SWIRE 
and UKIDSS-DEX (Deep Extragalactic survey), as well 
as partially covered by UKIDSS LAS; 1613 also covers a 
SWIRE (ELAIS-Nl) and UKIDSS-DEX field; 1645 cov- 
ers a SWIRE (ELAIS-N2) field. 2329 targets a DEEP2 
field. In addition to these targeting specific surveys, our 
survey fields have been used by other surveys. The Wig- 
gleZ survey (Drinkwater et al. 2010) has surveyed sev- 



TABLE 1 

RCS-2 PRIMARY SURVEY (PI IMAGING) 



Patch^ 



R.A. 



Dec' 



extent 



comments'^ 



RCS-2 primary survey (PI imaging) 



0047 


00:47.4 


4-00:45 


9x9 




0133 


01:33.2 


-00:10 


6x6 




0310 


03:10.3 


-14:11 


9x9 


Initial patcii extended to 










better cover the WMAP 










cold spot 


0357 


03:57.2 


-08:48 


6x6 




1040 


10:40.9 


-f-57:48 


6x6 


Lockman Hole 


1111 


11:12.0 


-05:52 


9x9 




1514 


15:14.6 


+08:55 


9x10 




1613 


16:13.6 


+55:00 


6x6 


ELAIS-Nl 


1645 


16:45.9 


+39:35 


6x6 


ELAIS-N2 


2143 


21:41.0 


-00:07 


10x10 




2329 


23:26.1 


-01:14 


8x6 


DEEP2 


2338 


23:38.5 


-09:07 


9x9 




CFHT Legacy Survey Wide component 


Wl 


02:18.0 


-07:27 


9x8 


XMM-LSS 


W2 


08:57.8 


-03:18 


5x5 




W3 


14:17.9 


+54:30 


6x7 




W4 


22:11.4 


+01:48 


6x6 


VVDS 



patch name 

RA & Dec (J2000) of the patch centre 

size (in number of 1 deg^ pointings) .This gives the approx- 
imate extent of the patch, but the layout is not necessarily 
rectangular. The exact geometry can be seen in Fig. 3. 

other relevant comments, such as overlapping surveys 

eral of our patches (0047, 0133, 2143, 2329, 2338) with 
GALEX in order to select z>0.5 Lyman Break Galaxy 
candidates using optical/UV colors to generate a large 
spectroscopic redshift sample to measure cosmological 
parameters via the signature of baryon acoustic oscilla- 
tions imprinted on their clustering. The Canada Prance 
High-z Quasar Survey (CFHQS, Willott et al. 2005) uses 
our z-band imaging in conjunction with their own z-band 
imaging to search for high redshift (2^7) quasars. 

These fields total 785 one-square degree pointings of 
imaging data for RCS-2. To this we add the four patches 
of the publicly- available Wide component of the CFHLS 
which tallies 171 cleg^ in u* , g' , r', i', z' . The reformat- 
ting and recalibration of these data to resemble the re- 
duction of our PI imaging data is described in §8. Thus, 
overall, RCS-2 comprises 955 one-square degree point- 
ings. 

3. PRE-PROCESSING 

Standard CCD pre-processing (bias subtraction, 
bad-pixel masking, flatfielding, etc.) is performed 
by the Elixir Project^ at CFHT. For the inter- 
ested reader, an excellent description of the raw 
MegaCam data, including example images at var- 
ious stages of Elixir reduction, may be found at 
http : //www . of ht . hawaii . edu/Science/CFHTLS-DATA/rawdata . h 

We use the default Elixir reduction to handle all these 
preliminary steps, except for the 2;-band defringing which 
we found could be improved upon. Due to the large vol- 
ume of data generated by our observations, rather than 
starting from raw (non-defringed) data and attempting 

® see http://www.cflit.hawaii.edu/Instruments/Elixir/ 
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rT 14 tl- 1.1 

Fig. 2. — Layout of a typical patch (0133). Points show objects in the photometric catalog brighter than r < 20. The chip gaps within 
each pointing are clearly visible. Pointings are labelled with their two character names, as described in the text. Overdensities of objects 
are visible where pointings overlap, as these duplicate areas have not yet been removed in the stitching process (see §9.2). Artefacts such as 
satellite trails have not yet been cleaned (see §9.1) and examples of these can be seen in pointings C5 and E3. Additional gaps are visible 
from where one of the CCDs (lower right corner) was non-functional when EO and FO were observed. 
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Fig. 3. — Layout of RCS-2 patches on the sky in Cartesian projection, plotted as a function of RA (hours) and Dec (degrees). Grey-scale 
shows distribution of Galactic Extinction from Schlegel et al. (1998). Squares denote individual Megacam pointings. (Horizontal gaps 
between pointings in high declination patches are an artefact of the projection and the pointings do in fact overlap, as described in the 
text.) Patches are labelled with patch names. Additional (unlabelled) pointings around 13 hours equatorial denote a patch which was left 
uncompleted (1303) and will not be considered as part of the survey proper (due to the non-contiguous nature of the observed regions), 
but individual pointings will still be included in searches for strong lenses. 
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to generate new master fringe frames from other z-band 
observations taken at similar times to our observations^, 
we investigated the possibihty of starting with the Elixir 
(non-optimally)-defringed z-band images. It was found 
that running these images through defringing code de- 
veloped by the CFHT Supernova Legacy Survey (SNLS) 
Team (kindly provided by A. Conley) resiilted in im- 
provement to the defringing for all z-band images.'^ How- 
ever, for a small fraction of the z-band data (^5% of all 
our images), significant fringing is still present. The ini- 
tial fringe amplitude in the raw data was at a level of 
'--^15% of sky*^. This was typically reduced to better than 
5% of sky by the initial Elixir defringing and reduced to 
negligible levels (^0.5% of sky) by the SNLS code for 
^95% of all our data. For the remaining ~5% of our z- 
band images, the residual fringing was at worst ~3% of 
sky. Thus, the most serious impact to our photometric 
measurements from the residual fringing would be an ad- 
ditional photometric error of f«0.03 mag in ~5% of the 
z-band photometry added in quadrature to the intrin- 
sic Poissonian error due to photon statistics. Thus, this 
inflates the photometric errors for very bright objects 
where the intrinsic error is ^0.03 mag; but has negligi- 
ble impact on the vast majority of our objects, which 
are fainter and thus have larger intrinsic errors. The 
one aspect of our catalog construction impacted by the 
presence of low-level residual fringing is object detection. 
The implications for this are discussed in §4. 

We note that our final images are identical to those 
generated by Elixir and stored in the CFHT archive (ex- 
cept for the 2;-band which have been further de- fringed). 
However, for each image, we also generate an associated 
bad pixel mask (BPM) using the procedures outlined be- 
low, which stores information about pixels contaminated 
with reflection haloes from bright stars, satellite trails, 
and cosmic ray hits, etc. 

3.1. Image alignment 

An initial estimate of the world coordinate system 
(WCS) is taken from the Elixir solution written in the 
image headers. In many cases this solution can be off- 
set from the real pointing position by several arcseconds, 
and a significant fraction of these show a non-physical in- 
strumental solution (such that one or more of the chips 
appear to be overlapping or displaced from their nominal 
position in the camera grid). So, we construct a corrected 
WCS by running fast object-finding, with a high thresh- 
old set to find only a subsample of the brightest objects, 
and pattern-match these against an astrometric reference 
catalog. We use a model for the global placement of the 
chips within the camera and solve for the locations of all 
objects within the pointing simultaneously. This allevi- 
ates the problem with displacement of individual chips, 
from the Elixir solution, which would otherwise cause 

'' due to the queue-scheduled nature of the observations, it is not 
always possible to find sufficient z-band data taken within a short 
enough time-frame to do this anyway. 

* Briefly, this code uses principal component analysis (PC A), 
taking a set of principal components (the eigenvectors) of the z- 
band fringe pattern from data taken on an early SNLS observing 
run (using McgaCam). It then attempts to find the amplitude of 
the fringe pattern wliieh must be subtracted (the eigenvalues) in 
order to remove the pattern. 

^ see http : //www . cf ht . hawaii . edu/Science/CFHTLS-DATA/rawdata 



problems when we attempt to measure offsets between 
the different filters for each pointing. This quick calibra- 
tion follows a similar method to that which provides our 
final, accurate astrometric calibration, which is described 
in detail in §6. For each pointing, a WCS is generated 
for each of the three filters, g, r, z, and this is used to 
find the approximate offset between the different filters. 

4. OBJECT DETECTION AND PHOTOMETRY 

Object finding and photometry arc performed using an 
updated version of Picture Processing Package (PPP; 
Yee 1991). While the basic methods are the same as 
those described in Yee (1991) and Yee et al. (1998), be- 
cause of the large amount of imaging data involved, a 
considerable number of additional algorithms have been 
implemented to allow the procedures to be carried out 
as a completely automated pipeline. In the following, 
we provide brief descriptions of the basic methodology 
and the added features that allow the automation of the 
process. 

4.1. Object finding 

In RCS-1, object detection was performed on coadded 
Rc + z images in order to increase the overall depth of 
the observations and to ensure the inclusion of the very 
reddest objects which may have been missed if the Kq- 
band alone (the deeper band) had been used. As a resiilt 
of the low-level (^3%) residual fringing in a small frac- 
tion (^5%) of the MegaCam z data, it is not possible to 
everywhere use co-added r -\- z images for object detec- 
tion. (The z-band fringing means that spurious objects 
are more likely to be detected if they sit near the peak of 
a fringe, whereas some real objects will be missed if they 
sit in the valley of the fringe pattern. The higher noise in 
the sky caused by fringing also artificially increases the 
detection threshold, meaning that some genuine objects 
may be missed in clean areas.) So, for consistency, and to 
be conservative, object detection is performed solely on 
the r-band images. Our depths in r and z are matched 
(by construction) for z~l red-sequence cluster galaxies, 
and so even though the catalog is formally r-selected it is 
also z-limited. Only for galaxies redder than t,^\ cluster 
members will the catalog begin to become incompleted^. 
Therefore, the only slight limitation is that searching for 
extremely red objects (r-dropouts) is not possible with 
this catalog. 

PPP performs object detection on each chip as de- 
scribed in Yee (1991), with additional algorithms to au- 
tomatically compute the detection threshold and to block 
out problem regions on the r-image. If artifacts (re- 
flection haloes and saturated columns) associated with 
bright stars contaminate the chip being considered, the 
measurement of the sky noise level, which is used in the 
estimate of the object detection level, may be biased 
high. Furthermore, these regions will also contain many 
false object detections. So, a mask is applied to remove 
contaminated areas from consideration when estimating 
the sky level and the noise. 

To mask bright stars and their halos, the approximate 
WCS solution derived in §3.1 (typically still accurate to 

this limit is tagged on a frame-by-frame basis by the PPP- 
. htrribasured photometric limits. 
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better than 1") is used to send a query to Vizier^'^ and 
download a catalog of all bright stars from the Tycho 
2 catalog (H0g et al. 2000) within the pointing. This 
bright star list is used to flag regions in which the sky 
is heavily contaminated by light from a bright star, and 
which must not be used when estimating the sky level 
for object detection. At this stage, a simple empirical 
relation between the observed size of the star halo and 
the cataloged magnitude of the star is used to set the 
mask size (which is set to be conservatively larger than 
might be needed). Later in the pipeline, a more detailed 
model of the properties of the reflection haloes around 
bright stars is used to delimit the area for object find- 
ing and to set mask flags in the final catalog (sec §9.1). 
The masking of saturated columns is done by identifying 
and connecting pixels above the saturation value starting 
from the positions of the bright stars. The masked areas 
are used neither in estimating detection threshold, nor 
in object finding. 

Object detection otherwise follows the same basic pro- 
cedure as described in Gladders & Yee (2005). Peaks 
are identified as significant enhancements measured in 
a 3x3 tapered box more than 2.6a above the local sky 
(excluding masked regions, as described above). Two cri- 
teria based on minimum mimbcr of connected pixels in a 
smoothed image and the 'sharpness' of the candidate in 
the unsmoothed image are used to eliminate cosmic ray 
and noise spike detections. 

4.2. Photometric measurement 

Total photometry of an object is derived based on the 
growth curve of the object which is measured in a se- 
ries of concentric circular apertures around the object, 
masking nearby objects as required (see Yee 1991 for de- 
tails). The growth curve is initially measured to 8.5" for 
all objects. The center of the apertures is determined 
by an iterative procedure that is accurate to fractional 
pixels (down to better than 1/10 of a pixel, depending 
on the signal-to-noise ratio of the detection). To allevi- 
ate the problem of pixelation at small radii, each pixel is 
subdivided into 7x7 subpixels before integrating the flux 
within an aperture. The photometric curve of growth so 
produced is used to identify an optimal aperture within 
which to measure the total magnitude following the pro- 
cedure described in Yee (1991). If the optimal aperture 
of the object is smaller than the standard aperture of 
8.5" diameter, the magnitude within the optimal aper- 
ture is extrapolated to the standard aperture using cor- 
rections derived from the shape of the growth curves of 
bright point-source objects. The standard aperture size 
is chosen to include close to 100% of the light for small 
objects under all our seeing conditions. For brighter re- 
solved objects, which normally would still have increas- 
ing flux at the 8.5" aperture, growth curves extended 
to a maximum diameter of 25" are measured and used 
in determining the optimal aperture to make sure that 
the bulk of the light is included. We note that objects 
with an optimal aperture smaller than 1.5" are classified 
as non-detections. The error on the total magnitude is 
then calculated as the sky noise within the photometric 
aperture used, which is sky noise-limited for the faint 
galaxies of interest here. 

see http://vizier.hia.nrc.ca/ 



The object positions detected from the r image are 
taken as the 'master' position list. To perform photom- 
etry on the images of the other filters, the position list 
has to be transformed to the pixel coordinates of these 
images to an accuracy of about 2 pixels, so that PPP 
is able to determine without ambiguity the subpixel cen- 
troids of the objects. For images that are offset from the 
parent r image by less than 10 pixels, or, images whose 
offset can be estimated to an accuracy of better than 
10 pixels from their approximate WCS solution (§3.1), a 
simple and very fast algorithm is used to register the po- 
sition file. Here, 35 to 70 brightest non-saturated objects 
are identified in the r image using a quick-find algorithm 
and used as position reference objects. Local maxima 
in the daughter frame within 10 pixels of the expected 
positions of the reference objects are located and used 
to derive average shifts in x and y. A rotation, if one 
exists and is sufficiently small (<~1 degree), is also mea- 
sured by comparing vector angles between pairs of ob- 
jects. This transformation is performed on the position 
list, and photometry carried out on the daughter frame. 

Occasionally, registration of the object positions for 
frames for the other filters of a pointing is not so straight- 
forward due to a significant rotation and even a small 
scale change. This could occur when MegaCam is re- 
moved and re-installcd, or the optics of the camera are 
changed (as in the case when one lens was inverted) or 
a filter is replaced (as in the case of the i-band filter, see 
§7). In these cases, more CPU-intensive algorithms are 
used to deal with the registration of the position files. 
The basic algorithm is a brute-force technique of find- 
ing the best X"Square fit of the two sets of bright po- 
sition reference objects, stepping through small incre- 
ments of xy shifts, rotation, and scale. The algorithm 
steps through increasingly more parameters in the above 
order, so as to minimize computational time. PPP per- 
forms the overall registration process by starting with 
the simplest procedure, tests the validity of the trans- 
form (based on whether the transformed positions from 
the r frame match the actual objects on the other filter 
frames close enough for PPP to compute an accurate 
fractional pixel centroid), and if not, performs the next 
more complicated matching procedure. In most cases, it 
is the i-band images (discussed in §7), sometimes taken 
a year or more apart from the others, and comprising 
images taken through two different filters over the whole 
survey, that require the more time-consuming registra- 
tion procedures. 

With growth curves for each object measured in all 
filters, colors are measured using a smaller aperture to 
improve the signal-to-noise ratio by minimizing the sky 
noise. The "color aperture" used is either a 3" diameter 
aperture or the optimal aperture, whichever is smaller. 
The same aperture is used for each filter, taking the r- 
band as the reference. No correction is made for seeing 
differences between the bands since the seeing values are 
typically comparable and the aperture is already rela- 
tively large compared to the PSF. The total magnitude 
for each of the other filters is then calculated by adding 
the relevant color to the total magnitude estimated in the 
reference filter (r) . The color error is then the quadrature 
sum of the errors in the two color apertures used to com- 
pute the color. Both color errors and total magnitude 
errors are propagated through to the final photometric 
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Fig. 4. — The distribution of seeing values (one value for each 
chip) in each filter. 

catalogs. 

A 5-cr reference limiting magnitude for each chip is es- 
timated by scaling the flux of a set of bright reference 
stars until they reach a signal-to-noise level of five. Typ- 
ically for extended sources, the 100% completeness limit 
is about 0.6 to 0.8 magnitudes brighter than this limit 
(Yee 1991). An estimate of the seeing is made by mea- 
suring the average FWHM of the reference PSF stars 
identified in each chip. The object catalogs typically 
comprise around 4000 objects per CCD chip. The val- 
ues of seeing measured by PPP, one value per chip, are 
shown in Fig. 4. The median values of the seeing are 
[g,r,i,z] = [0.79", 0.71", 0.53", 0.67"], but note that the 
distributions are broad due to our two-tier strategy. A 
significant fraction of our data possesses seeing better 
than ^ 0.6"in r. 

4.3. Star- Galaxy Classification 

Star-galaxy separation is performed by comparing the 
shape of the growth curve of each object to the weighted 
average growth curve from a set of four to eight refer- 
ence PSFs from nearby unsaturated stars (see Yee 1991 
and Yee et al. 1996) chosen from a list of automatically 
identified reference PSFs from the same chip. The list of 
identified reference PSFs is created from all bright unsat- 
urated objects that do not have close neighbors and are 
deemed 'not fuzzy' by an algorithm that iteratively elim- 
inates objects that are 'fuzzier' than the average. Each 
object is then given a classification of 0-artifact/spurious 
object ; 1 or 2-galaxy; 3-star; 4-saturated (Yee 1991). 
The subdivisions of types 1 and 2 for galaxies is for his- 
torical reasons related to the development of PPP and 
no distinction is made in practice. An object is classified 
as saturated if one or more pixels is above the saturation 
ADU value as provide in the Elixir FITS header of the 
image. 

Based on tests using simulated images, the star-galaxy 
classification is robust up to the 100% completeness limit 
under typical seeing conditions (Yee 1991). At magni- 
tudes fainter than 5a, star-galaxy classification is per- 
formed on a statistical basis using the "variable classi- 
fier" procedure based on the relative sharpness (rather 
than absolute sharpness) of the objects as described in 
Yee (1991). 

5. PHOTOMETRIC CALIBRATION 



Accurate photometric calibration is one of the most im- 
portant aspects of any imaging survey. For RCS-2, given 
the color-selected nature of our galaxy cluster sample, 
obtaining the highest possible accuracy in the calibra- 
tion of each color is particularly important. Here we 
describe the technique we employ to calibrate the colors 
of galaxies as accurately as possible using the colors of 
stars in our fields. 

The photometric catalogs generated by PPP all as- 
sume a single photometric zeropoint for all of the obser- 
vations made in each of the three filters. Some of our 
observations were made in mildly non-photometric con- 
ditions, a condition we allowed after initial tests with our 
color-calibration technique demonstrated that we could 
achieve accurate calibration without the use of standard 
stars. Only mildly non-photometric conditions were al- 
lowed so as to not unduly decrease the depth of the ob- 
servations, but this still increases the number of usable 
nights available to our project, since we do not require 
photometric conditions. 

In RCS-1, the [Rc — z') color calibration of each point- 
ing was performed by requiring the colors of bright galax- 
ies to agree with a reference pointing. The Rc-hand ze- 
ropoint was then set by requiring the number counts of 
galaxies to agree from pointing-to-pointing. In principal, 
the uniformity of the colors of stars or galaxies may be 
used. However, since we are primarily interested in ex- 
tragalactic studies with RCS-2, we find it preferable to 
force-fit the properties of stars to agree from pointing-to- 
pointing in our survey, decoupling any calibration errors 
from the properties of the galaxies we wish to study. For 
example, forcing the number counts of galaxies to agree 
from field to field then makes it impossible to study den- 
sity variations of galaxies on scales comparable to the size 
of each field, since we would have, by construction, forced 
all these values to agree. Thus in RCS-2 we will use the 
colors of stars to accurately set the calibration of each 
color, and then calibrate the value of the magnitudes in 
one reference band by comparison with an overlapping 
reference survey (2MASS). 

5.1. Fitting the stellar locus 

We use the uniformity of the colors of Galactic stars 
in order to set the color calibration from pointing-to- 
pointing in RCS-2. This method is an extension of the 
idea used by Hsieh et al. (2005) to calibrate multicolor 
{B, V) follow-up observations of RCS-1 fields. Hsieh 
et al. (2005) constructed histograms of objects classified 
as stellar and required that their {z' — Rc), (V — Rc) and 
(B — Rc) distributions agree from pointing-to-pointing 
(holding the i?c-band calibration fixed). 

Our method is inspired by Ivezic et al. (2004) who used 
the uniformity of the stellar locus in high-Galactic lati- 
tude fields in SDSS to assess the accuracy of their (inde- 
pendent) photometric calibration (using standard stars). 
They showed that, for high latitude (|&| > 15 degrees) 
fields, the colors of faint stars are essentially identical 
and arc dominated by stars in the Galactic halo (and 
thus behind the majority of the Milky Way dust). They 
showed that one can define principal colors along the stel- 
lar locus, and that constructing a histogram along these 
principal colors results in a Gaussian distribution where 
the width is dominated by the intrinsic width of the locus 
(~0.02 mag) if the photometric errors are small. Thus we 
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combine the idea of principal stellar colors for checking 
the accuracy of photometry (Ivczic ct al. 2004) with the 
color histogram-matching technique (Hsieh et al. 2005) 
to form a method for calibrating imaging surveys from 
the principal colors of stars. 

Our procedure is as follows. We begin by constructing 
a reference sample. Since our primary aim is a highly 
uniform color calibration across our entire survey, we de- 
sign a calibration set using only MegaCam data to cal- 
ibrate our filters internally to the native MegaCam sys- 
tem. We can then transform our (uniformly-calibrated) 
catalogs to any other (standard) system afterwards. To 
build a reference set, we select a subsample of around 20 
RCS-2 pointings which overlap with SDSS. Both datasets 
are corrected for Galactic Reddening using the map of 
(Schlegel et al. 1998), since the majority of the {\b\ > 15 
degrees) stars lie behind the Milky Way dust (Ivezic et al. 
2004). (All our fields lie at |6| > 30 degrees, with the ex- 
ception of the CFHT Legacy W2 field which is |&| > 20 
degrees.) For each pointing wc match individual objects 
with their counterparts in SDSS and work out the me- 
dian offset from SDSS for that pointing and filter. For 
example, for the i*'* pointing in the r-band filter, the 
natively calibrated magnitude, ri^cAL, is related to the 
instrumental magnitude produced by PPP, ri^instr, by 
the offset dvi by 



where 



fi,CAL — fi,instr "I" 



dn = median{rij^sDSS ~ njjnst 



(1) 



(2) 



and Vij^sDSS and rij^mst refer to the j*" object matched 
between SDSS and MegaCam respectively. Now, for the 
r-band, the color term between MegaCam r and SDSS r 
is negligible, so this transformation effectively transforms 
onto the SDSS AB system. For the other filters there 
are color terms. Adopting the above procedure for all 
filters (i.e., neglecting to specifically apply a color term) 
means that our reference set is calibrated onto a native 
MegaCam system which agrees with SDSS for objects 
with the colors of the median-color objects in SDSS. We 
derive the color terms which may be applied to convert 
the RCS-2 native system to SDSS AB at the end of this 
section. 

We can now take each (uncalibrated) science point- 
ing in turn and calibrate it to agree with the above set. 
We do this for every pointing, including those which were 
used in constructing the reference sample, to ensure com- 
plete consistency across the survey. For calibrating g, r, 
z data we will calibrate the {g — r) and (r — z) colors, 
initially holding the r-band magnitude fixed. We start 
with the (g — r) color since the red branch of the stellar 
locus is almost independent of (r — z) color (Fig. 5, see 
also fig. 2 of Ivezic et al. 2004). 

The {g — r) principal color is obtained by simply sub- 
tracting the {g — r) color of the peak of the [g — r) his- 
torgam, which we denote {g~r)o, where {g — r)o = 1.289. 
The resulting histogram, now centred on zero color, is 
shown as the upper right panel in Fig. 5. 

The principal color used to calibrate (r — z) is defined 
by fitting a linear relation to the blue part of the color- 
color diagram (Fig. 5): 

PI = {r-z)-[0.779{g-r)-0.152]n{g-r) < 1.158. (3) 
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Fig. 5. — Main panel: color-color diagram of a subsample of stars 
used in the reference dataset. Dashed lines indicate the positions 
of the principal colors used. The vertical line is the (g — r) principal 
color which is simply (g — r) — (g — r'o) where (g — r)() is the reference 
color. A histogram of this principal color is shown in the upper 
right panel. The lower right panel shows a histogram of the PI 
principal color, described in the text and indicated as the sloping 
dashed line in the main panel. See text for details. 

The (g — r) color cut ensures that only the blue branch 
of the stellar locus is used. This results in a very clean 
Gaussian in the histogram, shown as the lower right panel 
in Fig. 5. 

5.2. Application of the fitting procedure 

With the above reference set in-hand, each science 
pointing may be examined in turn and its star colors 
forced to agree with that of our reference set. All ob- 
jects classified as unsaturated stars with magnitude er- 
rors < 0.2 mag are extracted and color histograms built 
in the same way as for the reference sample. The {g — r) 
histogram in the science pointing is cross-correlated with 
the of the reference pointing, and the offset A{g — r)i 
found. This procedure is typically found to be accurate 
to a few hundredths of a magnitude (w 0.05), by com- 
paring the result with an independent calibration such as 
SDSS. In order to improve this accuracy further, a Gaus- 
sian profile is fitted iteratively around the position of the 
peak of the histogram, and the peak position and width 
recorded. The offset between the Gaussian-fitted peak 
and zero is a refined estimate of the color offset, A(g— r)2. 



The width of the fitted Gaussian, a. 



ia-r) 



is used to check 



for problems with the photometry/calibration method. 
The width of the distribution is set by the intrinsic width 
of the stellar locus convolved with our photometric er- 
rors, and thus should be approximately constant from 
pointing-to-pointing. If iy{g-r) > 0.15 mag then the 
pointing is flagged as having potential problems with the 
photometry. The inflated width of the Gaussian usu- 
ally means that the {g — r) colors are smeared out due 
to the the g and r frames not having been registered 
correctly (occasionally due to the wrong pointing posi- 
tion being observed by the queue observers). A similar 
warning flag can be generated by flagging unusually low 
amplitudes of the Gaussian peak (since the number of 
stars in each pointing is approximately constant). Such 
calibration problems are later inspected manually. 
The calibration in {g — r) is thus 

(g - r)cAL = (5 - r)^nst + A(5 - r)i + A{g - r)2, (4) 
where A{g — r)i is the offset determined by cross- 



10 



Gilbank et al. 



correlation of the color histograms, and A(g — r)2 is 
the offset from the peak of the itcratively-fitted Gaus- 
sian. In the case of a Gaussian width cr(g-r) > 0.15 mag, 
A{g — r)2 is set to zero and a warning flag set. In a small 
minority of cases, this warning flag is caused by some 
problem other than bad registration of the images (such 
as contamination by satellite trails). In this case, if vi- 
sual inspection of the calibration histograms reveals no 
obvious problem, then the calibration is accepted using 
just the cross-correlation offset, A{g — r)i. 

A similar procedure is applied in (r — z), but using the 
principal color PI (Eqn. 3) when constructing the color 
histogram. Note that the {g — r) color cut used with PI 
to select the blue branch of the stars gives a very sharp, 
clean histogram in {r — z). This can be done accurately 
since the {g — r) color has been calibrated in the pre- 
vious step. Again, a cross-correlation offset followed by 
a Gaussian fit to the located peak is performed and the 
same check, an < 0.15 mag, made of the results. 

5.3. Magnitude calibration 

The above procedure results in very accurate calibra- 
tion of the {g — r) and (r — z) colors (^0.03 mag). How- 
ever, the magnitude in any individual filter has not yet 
been accurately calibrated and may be considerably in 
error by several tenths of a magnitude (the data may 
have been observed through thin cloud, and even if not, 
we have made no explicit adjustment for atmospheric ex- 
tinction). Operationally we held the r-band fixed above, 
but in order to calibrate the individual magnitudes it is 
sufficient to calibrate any one of the individual filters, 
provided we hold the colors fixed to their calibrated val- 
ues. To obtain a calibration, we need a survey which 
overlaps with objects in our sample and possesses well- 
calibrated photometry. SDSS would be ideal, but it does 
not overlap with all of our pointings (or even all of our 
patches). So, we use NIR data from 2MASSi2. J-band 
total magnitudes from the 2MASS point source catalog 
are downloaded for objects in the RCS-2 footprint. A 
reference set is created again using the SDSS-calibrated 
pointings, as above. Thus, we can construct color-color 
diagrams expected for stars in any combination of our 
(-I-2MASS) observed filters. We use (g — J) versus (g — r) 
as shown in Fig. 6. As mentioned above, we only need to 
calibrate one filter. Although z-band is closer in wave- 
length to observed J-band, it does not matter which com- 
bination is used, since the fit is performed in color-color 
space and the locus is just as well defined in these col- 
ors. 

For (g — J) the principal color used is: 

P2 = {g~J)- [2.226(.g - r) - 0.764] n [g - r) < 1.2. (5) 

The offsets, calculated as above, result in a g-band 
shift, Ag, which must be applied to each filter in such a 
way as to preserve the calibrated colors; i.e., if we re- write 
the original color shift (Eqn. 4) as shifts to the individual 

We thank our WiggleZ collaborators for originally suggesting 
this idea. 

The reason for choosing g is that this is the last of the three 
filters on which PPP performs photometry, so using it as the cali- 
bration filter is useful for catching problems with the photometry 
— if the last filter (g) fails, then all the photometry fails for that 
pointing. 
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Fig. 6. The reference pointing used to measure the (g — J) — 
— {g — r) locus by combining 2MASS and RCS-2 photometry. The 
blue branch in [g — r) exhibits the smallest scatter and so this is 
fitted with the function in Eqn. 5 (dashed line). 

filter, defining the r-band shift as zero, Ag = A{g — 
r) and Ar — 0. The calibration offset measured from 
2MASS, Ag2MASSj gives resulting overall corrections to 
the instrumental magnitudes of 



and 



Afiftot = Ag + Ag2MASS 



Artot = Ar + Ag2MASS 



(6) 



(7) 



to preserve the calibrated {g — r) color, and where Ar 
was above defined as zero. 

Similarly, for z the offset would be 



AzTot = Az + Ag2MASSi 



(8) 



since its calibration was performed in (r — z) and r has 
just had Ag2MASS added to it. 

Just to reiterate the procedure, we take data, which 
may have been acquired in non-photometric conditions, 
and obtain an accurate photometric calibration simply 
by forcing the colors of stars to agree with each other 
from pointing-to-pointing (which ensures highly uniform 
color calibration) and obtain a magnitude calibration by 
forcing these colors to agree with a reference set including 
J-band magnitudes from 2MASS in one of the colors. 
This procedure negates any need for separate correction 
of observational effects such as exposure time differences 
and atmospheric extinction. We note that the idea of 
using stars to provide or test a photometric calibration is 
not a new one (indeed. High et al. 2009 recently presented 
an implementation of a similar technique) , but we believe 
we are the first to implement such a method as the sole 
means of calibrating a survey of this size. 

It is worth emphasising that since the stars used in the 
stellar locus-fitting arc behind the majority of the Galac- 
tic dust, and we have reddening-corrected our reference 
dataset, the magnitude offsets measured between instru- 
mental magnitudes and the reference set will correct the 
data to a system already corrected for Galactic Extinc- 
tion. We note that in practice we prefer to both redden- 
ing correct our instrumental science catalogs and our ref- 
erence dataset before performing the calibration. Thus, 
the magnitude differences, Amag (e.g., Eqn. 4), can be 
directly applied to the initial instrumental magnitude 
catalogs (i.e., not corrected for Galactic Extinction) and 
this corrects the catalog initially to observed rather than 
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Fig. 7. — 5(7 point source limiting magnitudes in each passband 
as measured by PPP. 

reddening-corrected magnitudes. We find this more use- 
ful for considering observational effects (such as checking 
image depth). We store sub-pixel interpolated versions 
of the Schlegel et al. (1998) dust maps (to prevent pos- 
sible sharp discontinuities across our patches caused by 
the large pixel size of the former) matched to our survey 
regions as part of this process. These then allow straight- 
forward reddening correction of the observed magnitudes 
afterwards. 

A second useful note about our procedure is that it is 
the almost vertical nature of the red branch of the stellar 
locus in {g — r) — {r — z) which greatly simplifies the initial 
calibration. 

The final depths of our imaging data are shown in 
Fig. 7. Limits refer to the 5cr point source limiting mag- 
nitude as measured by PPP. 

5.4. Comparison of calibrated data with SDSS 

In order to assess the accuracy of RCS-2 photome- 
try, the significant overlap with the large, well-calibrated 
SDSS is used. In addition, this allows calculation of the 
transformation between the RCS2 MegaCam native sys- 
tem and the widely used SDSS system. 

The first step is to perform object-by-object matching 
for each pointing, as was done to construct the reference 
stellar locus dataset, as described in the previous section. 
To elucidate this, we show the comparison of a randomly 
selected RCS-2 pointing in Fig. 8. This shows the com- 
parison between r-band magnitudes in RCS-2 and SDSS 
for this pointing. The average offset between the two sets 
of calibrated magnitudes is estimated via a biweight fit 
to the difference which is shown by the solid horizontal 
line. For this example pointing, the difference is 0.048 
mag. The scatter is dominated by SDSS photometric er- 
rors at the faint end, due to the shallower SDSS photom- 
etry. At the bright end, the scatter is likely dominated 
by differences in the way magnitudes are measured by 
each survey. As described above, the RCS-2 magnitudes 
from PPP are essentially estimates of the total magni- 
tude, and these are compared with SDSS ModelMag mag- 
nitudes. 

Comparing colors relies on a more stable measurement 
than trying to estimate total magnitudes, so we begin 
here. Fig. 9 shows the differences in (g — r) and (r — z) 
and (r — i)^"* colors measured in RCS-2 compared with 

the i-band data are described in §7. 



Fig. 8. — Object-by-object comparison of r-band total magni- 
tudes from a single, randomly-selected RCS-2 pointing with ob- 
jects in common with SDSS. Only objects unsaturated in RCS-2 
are considered. The solid line shows the best fit to the >3600 ob- 
jects in common. The dashed line shows zero offset for comparison. 
The offset is dominated by the calibration uncertainty in RCS-2, 
assuming a calibration error of AmsiO.Ol mag in SDSS (Ivezic et al. 
2004). The scatter is dominated by photometric uncertainties in 
SDSS at the faint end, due to the shallower depth of SDSS relative 
to RCS-2. 

SDSS. g-band magnitudes have been corrected for the 
color term between MegaCam and SDSS g, as will be de- 
scribed below. Each panel shows a different RCS-2 patch 
and each value in the histogram represents the biweight 
average offset (as discussed above) difference between the 
two surveys for a single pointing. Only pointings with a 
significant number of objects (typically >1000) in com- 
mon are considered. For each patch, the histogram in 
a given color is relatively narrow (;^0.05 mag), implying 
that the color calibration produces very uniform colors. 
There appear to be some small systematic offsets with 
respect to the SDSS photometry in the sense that the 
{g — r) colors are systematically marginally redder than 
SDSS (~0.03 mag), whereas the [r — i) and {r — z) colors 
are slightly bluer (~0.04 mag). There also appears to 
be some evidence that the size of this small offset varies 
from patch to patch. This might be due to inaccuracies 
in the assumption that the Milky Way stellar populations 
used in the color calibration are identical for every patch 
(either due to intrinsic differences due to, e.g., metaUic- 
ity of the stars, or inaccuracies in the assumption that 
they may all be considered behind all the Milky Way 
dust, as measured by the Schlegel et al. 1998 dust map). 
Nevertheless, the calibration accuracy of the colors still 
represents a significant improvement over that usually 
achievable with the typical number of standard star ob- 
servations made throughout an observing run. This is 
even more remarkable considering many of our data were 
taken in slightly non-photometric conditions. 

Fig. 10 is constructed in the same way as Fig. 9, but 
the histograms now show the difference in r-band magni- 
tude between RCS-2 and SDSS. The scatter is somewhat 
larger than the error in the color calibration, more like 
'-~^0.05 mag rms, with some larger outliers. This likely 
reflects the limitation of using 2MASS photometry to set 
the magnitude zeropoint. It is difficult to find a sufficient 
number of common stars between 2MASS and typical 4- 
m telescope photometry, such that the stars are unsat- 
urated in the larger telescope and sufficiently bright in 
2MASS to yield a high S/N measurement of the star's 
magnitude. Even with the one square degree field of 
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Fig. 9. — Check of the color differences from pointing-to-pointing by comparing RCS2 colors with those of SDSS. Amag = coliicS2 — 
colgoss- Differences are measured object-by-object (as in Fig. 8) and then a mean offset is assigned to each pointing. These histograms 
show the mean value of the offset for each pointing with significant SDSS overlap. Line styles denote different colors as annotated: (g — r) 
dashed blue line, (r — i) dotted red line, (r — z) solid purple line. 



MegaCam, only «30 stars are typically suitable for use 
in the magnitude calibration. 

5.4.1. Colour terms with respect to the SDSS system 

Differences in the shape of the response of the Mega- 
Cam system (filter -|-CCD-|-optics+atmosphere) versus 
SDSS can be readily calculated following an extension of 
the above comparison. Individually-matched objects are 
again considered, but this time the average offset in each 
filter for each pointing is first subtracted (i.e., this would 
make the histograms in Fig. 9 and 10 all zero- valued) . 
We do this because we are only interested in the shape of 
the response, i.e., the slope of the magnitude difference 
as a function of MegaCam color, not its absolute value. 
Fig. 11 shows this difference as a function of MegaCam 
{g — r) color. All common, unsaturated objects (>40 000) 
in a single patch are used. As can clearly be seen, the 
g-hand magnitude difference between the two surveys is 
a strong function of {g — r) color. The best fit relation, 
shown as the solid line, is given by 

5SDSS= 5 + 0.23 - 0.174(5 -r). (9) 

As can be seen from Fig. 11, the best- fit line passes 
through zero offset from SDSS at the position of the 
median color [{g — r) ^ 1.6], by construction, from 
Eqn. 2. Due to the way in which we have internally 
color-calibrated RCS-2, this transformation to the SDSS 
system is now a trivial last step of the calibration pro- 
cedure, which may be applied by using Eqn. 9 on the 
photometric catalog, if desired. 



Tests for color terms in r and z are performed in the 
same way. The results are shown in Fig. 12. The z-band 
shows a tiny residual color dependence « 0.01(r — z) 
which would make <0.01 mag difference to the z-band 
data over the typical color range of objects and so may 
safely be neglected when comparing with SDSS photom- 
etry. The r-band shows no measurable color term. 

5.5. Summary of photometric accuracy 

To summarise the systematic errors of our photomet- 
ric catalogues, our color calibration technique is able to 
produce colors calibrated to to an absolute accuracy of 
<^ 0.03 mag. This is performed independently of the 
magnitude calibration (single filter) which is somewhat 
less accurate (<; 0.10 mag) due to the need to tie to the 
2MASS system, which has limited numbers of objects in 
each of our pointings. The 5-a point source depths of 
our catalogs (at which the magnitude error is «0.2 mag) 
are given in Fig. 7. Our random magnitude errors are 
then limited by photon statistics in the usual way. One 
small exception to this is due to a small fraction (^5%) 
of our z-band data which possesses an additional error 
due to residual fringing at a level of <;0.03 mag. This 
makes a negligible additional contribution to the photo- 
metric errors for the vast majority of sources of interest 
to us for cluster detection (which are predominantly faint 
and have larger intrinsic photometric errors than this al- 
ready) . 

As one final test of the internal accuracy of our cata- 
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Fig. 10. — As for Fig. 9 but histograms now show the mean r-band magnitude offset for each pointing with respect to SDSS. 




Fig. 11. — g color term in {g — r) deduced from a comparison 
with SDSS g photometry. Plot shows MegaCam g filter compared 
with SDSS, g — gsDSSt ^ function of MegaCam (g — r) color for 
a whole RCS-2 patch. Solid line is color term equation given by 
Eqn. 9. 

logs, we present a comparison of photometric differences 

between objc^cts measured in the overlapping regions of 
neighboring pointings. Fig. 13 shows the results for one 
representative patch. 2143. For each pointing, the differ- 
ence between magnitudes in each of the g, r, and z filters 
measured within that pointing and one of its neighbors 
is recorded. Only pairs of pointings with at least 100 ob- 
jects in common are considered. The median difference 
between this pair of neighbors is then calculated and ap- 
pears as a single entry in the histogram in the left panel 
of Fig. 13. The median absolute deviation is recorded 
for the same set of objects and is entered into the his- 
togram in the right panel. The results are plotted for 



every unique pair of overlaps. The rms width of the dif- 
ference histogram (left panel) is better than 0.04 mag, in 
good agreement with the estimates from the SDSS com- 
parison. Similarly, the average intrinsic scatter within 
the overlap regions agrees with this. Tails toward greater 
values of the magnitude difference are seen, and illustrate 
one of the drawbacks of this comparison. The overlap- 
ping regions must necessarily come from the edges of 
the camera, where corrections for flatfielding, etc. are 
generally worst. Thus, comparing the edge in one point- 
ing with a different edge of the camera in a neighboring 
pointing generally produces the worst agreement, and the 
internal accuracy of our catalogs are likely better than 
this over most of the field. Nevertheless, this is a use- 
ful, if somewhat pessimistic, estimate of the uncertainty. 
Similar results are seen for the other patches. 

6. ASTROMETRIC CALIBRATION 

Astrometry is performed in a manner similar to that 
used in Gladders & Yee (2005). An instrumental solution 

is first generated using an observation of an astrometric 
standard field, taken as part of the MegaCam calibra- 
tion procedure. This is used to map individual chip co- 
ordinates to a global instrumental solution in terms of 
ARA, ADec from the pointing centre. This is achieved 
semi-interactively using custom-written routines to fit a 
2D cubic function using the IDL's POLYWARP routine. 
In the absence of differential atmospheric refraction and 
other real-world limitations such as the reproducibility 
of mounting the camera on the telescope, this solution 
would be sufficient to describe the relative positions in 
every observation. In practice, it is necessary to fit for 
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Fig. 12. — Estimates of r-band and z-band color terms with respect to the SDSS system from a comparison with SDSS photometry as 
a function of (r — 2) color for a whole RCS-2 patch. The solid line shows the best fit. For the z-band, the color term is 0.01 (which we 
ignore) and the r-band term is negligible. 
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Fig. 13. — Internal accuracy of photometry estimated from all overlapping pointings within one patch (2143). The left panel shows 
the median difference for each pointing compared with its neighbor, and the right panel shows the median absolute deviation of these 
overlapping objects. Each entry in the histogram is the average of all points in the overlap region between a unique pair of pointings. Only 
overlaps containing more than 100 objects are considered. See text for discussion. 



not only relative shift (of the pointing centre relative to 
the expected telescope pointing position) and rotation, 
but also for the differential displacement of objects across 
the field. This is achieved using lower-order corrections 
to the instrumental polynomial map, again using IDL's 

POLYWARP. 

For each pointing, x,y positions are taken from the 
PPP catalog. Now that the photometric calibration 
has been performed, it is a simple matter to extract a 
magnitude-selected sample in a similar way to an exter- 
nal astrometric reference catalog. The nominal object 
selection isl7<r<21, which is matched against the 
full depth of the USNO-Bl astrometric catalog (Monet 
1998) using the instrumental solution as an initial esti- 
mate of the positions. Shifts, rotation and distortion 
are found iteratively by proceeding in small steps. After 
each iteration, the rms offset from the USNO-Bl is mea- 
sured and the results recorded after the final iteration. 
A comparison of object positions from an entire RCS- 

With a survey as large as RCS-2, even relatively rare occur- 
rences can lead to a significant number of failures in automated 
routines. So, in practice we find it necessary to use a range of 
magnitude cuts if the nominal values fail to yield a sensible astro- 
metric solution. In very rare cases (the presence of a rich globular 
cluster, M2, and a nearby dwarf irregular galaxy, IC1613) it was 
necessary to intervene and reset parameters manually. 



2 patch with positions determined from SDSS is shown 
in Fig. 14. A generous matching radius of 3" is used be- 
tween the two catalogs to ensure that the difference is not 
underestimated due to distant matching objects being 
missed. The distribution of differences between object 
centres is shown with a Gaussian of mean, /i = 0.30"and 
width, a = 0.08". The distribution has a higher mean 
offset than might be expected given its width (the mean 
and variance should be the same if the distribution is 
approximately Poissonian) , and so the bulk of the offset 
may be due to systematic differences between the SDSS 
system and ours. A comparison between independent 
observations of the same objects in overlapping regions 
at the edges of our pointings shows that the rms differ- 
ence is closer to 0.1". Regardless, our astrometric pre- 
cision is ^ 0.3", and this is comparable to (or better 
than) the accuracy of the photographic plates used in 
making the USNO-Bl. This accuracy is confirmed by a 
large amount of spectroscopic follow-up with the wide- 
field spectrograph AAOmega by the WiggleZ collabora- 
tion (Drinkwater et al. 2010), and by our own follow-up 
observations. 

7. ADDITION OF i-BAND DATA 
7.1. Image registration and photometry 
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Fig. 14. — Astrometric offset between positions of SDSS objects 
and RCS-2 positions for all overlapping objects in a whole RCS-2 
patch (2329) . There are almost 60 000 galaxies in common. Dashed 
curve shows a Gaussian of = 0.30" and u = 0.08". 

Although the RCS-2 survey proper comprises g, r, z 
imaging, the majority of fields now also have i-band 
imaging obtained for the CFHQS (Willott et al. 2005) 

by reobserving our fields in i. We now discuss how these 
data have been incorporated, where available. The basic 
procedure is similar to that of the other filters, registering 
each z-band image to the completed grz catalogs using 
the alignment method described in §3.1. Once registered, 
photometry could proceed through ppp as before. 

7.2. Photometric calibration 

Due to the breakage of the original MegaCam i-band 
filter (i.MP9701) in June 2007, the survey contains i- 
band imaging from two different filters. The original fil- 
ter was replaced with a new one (i.MP9702) possessing 
a much bluer response (close to the SDSS response). As 
a result, a significant color term exists between the two 
filters. For simplicity in the photometric catalogs, we de- 
cided to convert the i.MP9701 native photometry to that 
of the i.MP9702 native response. This has the advantage 
that the latter has a negligible color term compared with 
SDSS i, and thus our i-band data will be automatically 
calibrated to the AB system. The transformation be- 
tween the two filters is found by comparing each with 
overlapping SDSS photometry using pointings in com- 
mon, as described for the other filters in §5.4. This is 
shown in Fig. 15. 

The conversion is 



= ioid + 0-0Q7{r-ioid), 



(10) 



filters respectively. The rms scatter in this transforma- 
tion is 0.034 mag. 

Colour calibration of individual science pointings is 
achieved in the same way as for the 2;-band. Stars with 
magnitude errors <0.2 mag are selected from each point- 
ing, and the stellar locus fitted to a reference locus in 
(g — r)-(r — «) , constructed as for the (g — r) (r — z) case. 
The principal color used (plotted in Fig. 16) is 

= (r-i)-[0.529(5-r)-0.125]n(g-r) < 1.158 (11) 

Since the other bands have already been calibrated, a 
single offset to the i-band is all that is required. 

8. CFHTLS-WIDE DATA 




Fig. 15. — ioid color term in (r — i) deduced from a comparison 
with SDSS i photometry. Plot shows MegaCam i.MP9701 filter 
compared with SDSS, i^id — isDSS: ^■s a function of MegaCam 
(r — i) color for a whole RCS-2 patch. Solid line is color term 
equation given by Eqn. 10. 




Fig. 16. — Colour-color diagram of the subsample of stars used 

as the reference set in the i-band calibration. Dashed lines indicate 
the principal colors used. The inset shows the histogram about the 
principal relation Pi (Eqn. 11), sloping dashed line. 

The CFHT Legacy Survey^^ is a joint project between 
Canada and France to image a number of fields using 
MegaCam. The 'Wide' component of the survey, con- 
sidered here, covers 171 square degrees in ugriz. The 
observing strategy is somewhat different from RCS-2, 
since multiple exposures with dithering arc used at each 
pointing position. Catalogues available to the Canadian 
and French communities (and later to the world) are 
produced by TERAPIX (Traitement Elementaire, Re- 
duction et Analyse des PIXcls de megacam)'^^. We use 
the T0006 internal release (2009 November) of the data. 
Briefiy, the TERAPIX reduction involves similar pre- 
processing to the Elixir products, with the additional 
steps of determining astrometric solutions and combining 
the dithered images in each filter in each pointing into a 
stacked image. Then object detection and photometry is 
performed using SExtractor (Bertin & Arnouts 1996). 

In order to make the TERAPIX catalogs resemble our 
method as closely as possible, we extract total magni- 
tudes in the r-band using SExTRACTOR's MAG_AUTO 
estimate. We then construct total magnitudes in the g, 
i and z bands by constructing the relevant color with 
respect to r-band using aperture magnitudes measured 

see http://www.cfht.hawaii.6du/Science/CFHLS/ 
see http : //terapix . iap . f r 
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in 4.4"diameter apertures and add these to the r-band 
total magnitude to mimic the method used by PPP. i.e., 

9Tot = {9- r)ap + rxot (12) 

for the g-band. where the subscripts Tot and ap de- 
note total and aperture magnitudes, measured as de- 
scribed, respectively. Star/galaxy classification is per- 
formed using SExtractor's CLASS_STAR parameter. 
TERAPIX flag all objects brighter than i < 17.5 as satu- 
rated regardless of whether they are actually saturated or 
not. To recover bright galaxies, which arc not saturated, 
wc follow the method of Lu et al. (2009) considering the 
half-light radius estimated by S Extractor as a func- 
tion of irot- 

The most important step in ensuring that the CFHTLS 
data and RCS-2 are on a uniform system is to ensure that 
the color calibrations arc the same, i-band data are first 
converted from the old to the new filter system using the 
color term described in §7.2. 

The above procedures produce catalogs in which each 
pointing is independently calibrated in terms of photom- 
etry and astromctry. The final steps arc to clean the cat- 
alogs, removing artefacts such as refiection halos of bright 
stars, satellite and meteor trails, etc., and to stitch the 
individual pointings in a patch into a contiguous area, re- 
moving duplicate objects in the overlaps between point- 
ings (as described in §9.2). The g, r, i, z data are then 
run through the RCS-2 photometric calibration pipeline, 
as described in §§5, 7.2. Since the TERAPIX catalog 
construction already contains a photometric calibration, 
this typically results in small ^ 0.05 mag shifts in the ze- 
ropoint of each filter in each pointing. TERAPIX provide 
two alternate photometric calibrations for their catalogs. 
The default involves calibrating a single pointing in each 
of their patches and then using this as a reference. The 
calibration of all other pointings in the patch is then de- 
termined via the overlaps of all overlapping pointings, 
holding the single reference fixed. TERAPIX also pro- 
vide photometric offsets which have been determined via 
Stellar Locus Regression using a technique similar, but 
not identical, to ours. In order to ensure complete unifor- 
mity across RCS-2, we use our stellar locus fitting and 
2MASS comparison to calibrate the colors and magni- 
tudes respectively, completely independently of the Elixir 
calibration. We find from comparison of each calibra- 
tion with individual objects in SDSS (as in §5.4) that 
our color calibration is somewhat more accurate than 
the TERAPIX color calibration (w0.02 mag vs «0.05 
mag) but that our magnitude calibration is marginally 
worse («0.05 mag vs «0.03 mag), again likely limited 
by the number of stars available from the overlap with 
2MASS and their associated photometric errors (see Lu 
et al., in prep for a detailed comparison). It is possible 
that we will modify the final RCS-2 magnitude calibra- 
tion to use internal overlaps as a final correction, but the 
current calibration is certainly sufficiently precise for our 
purposes. 

9. CLEAN CATALOG CONSTRUCTION 

9.1. Artefact masking 

Two distinct sets of artefacts are present in the data 
which must be masked out to produce a clean photo- 
metric catalog. The first set is predictable artefacts such 



as those associated with bright stars. The positions of 
bright stars in MegaCam pointings are known a priori 
(once the astrometric calibration has been performed), 
and so the locations of halo reflection artefacts may be 
automatically predicted and masked. Secondary reflec- 
tions within the camera lenses, or returned reflections off 
the detector or camera dewar window and then off an air- 
glass interface in the camera result in extended reflection 
haloes visible around bright stars. Careful examination 
of the images revealed three dominant reflection haloes. 
Each halo is of similar total intensity, and so the largest of 
these is visible above sky noise for only the brightest stars 
(typically at most a few stars per pointing) and the small- 
est is visible around relatively fainter stars (typically up 
to a few stars per mosaic chip). The outer boundaries of 
these artefacts are approximately circular, but the center 
of the reflection haloes of each type is a function of posi- 
tion within the camera reference frame. The halo sizes, 
and positions as a function of camera coordinates, were 
determined empirically from a set of test observations in 
each band, and these data were used to construct masks 
around known bright stars. The masking additionally 
requires a threshold brightness for each reflection halo; 
this was again determined empirically by examining ob- 
ject counts radially around bright stars, and choosing 
limiting magnitudes for each reflection halo that ensured 
that regions which are significantly compromised (either 
by having too many objects detected due to structure in 
the reflection halo, or by having too few objects due to 
the increased sky noise) are masked. Objects in these 
masked regions are retained in the photometric catalog, 
but a flag is set warning that they are located in a re- 
flection artefact and thus they may be straightforwardly 
rejected from subsequent analysis. 

The second set of artefacts are those whose locations 
cannot be predicted prior to examining the data, such 
as satellite trails and meteor trails (see Fig. 2). These 
are dealt with in a semi-interactive manner. Since these 
are linear artefacts, line-flnding techniques akin to the 
Hough transform (converting the positions of objects in 
the PPP catalog to polar coordinates and looking for 
overdensities of points in this coordinate space) can be 
used to flag potential artefacts. In practice, the images 
are examined interactively with sites of potential arte- 
facts flagged by this technique. At this point other arte- 
facts which may have passed the automated masking are 
also flagged. These include occasional reflection haloes 
around stars whose as-observed magnitudes are signifl- 
cantly brighter than their cataloged values. Also, bright 
and very low redshift galaxies - typically subtending tens 
of arcseconds or more - will tend to produce many ob- 
jects due to substructure in the galaxy. These regions 
are also masked to mitigate any influence of such objects 
on future analyses of object counts or clustering. Finally 
the most common masking added by hand at this stage 
is of diffraction spikes which have escaped the automated 
masking done by PPP (see §4). Typically these are ei- 
ther from close but resolved double stars which appear 
as single entries in the catalogs used to initially mask the 
images and which produce two sets of diffraction spikes 
neither of which is caught by the initial mask, or diffrac- 
tion spikes from stars which fall in the mosaic chip gaps. 
After all additional masking is complete, a mask image is 
constructed for each chip, with the values of the masked 
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pixels set to indicate the source of the masking. These 
images also include a record of missing filter data along 
image edges (due to misalignments between images of 
differing filters) with missing data in each filter flagged 
by a different mask value. In principle, this allows one 
to select subsets of data masked only for certain reasons, 
while including other nominally masked objects; this may 
be of interest under some circumstances. 

9.2. Patch stitching and random catalogs 

Up until this point, the data processing has been per- 
formed on individual pointings. Due to the survey strat- 
egy, overlaps occur between neighboring pointings. The 
final step is to 'stitch' the overlapping pointings into sin- 
gle contiguous patches, removing duplicate objects which 
occur in the overlapping regions. This is performed in 
the same manner as for RCS-1. Briefly, the midpoint 
of the joint area between neighboring pointings is found 
and then each field is truncated at the midpoint. See 
Gladders & Yee (2005) for details. Thus, only objects in 
unique, non-overlapping areas are retained in the stitched 
photometric catalog. 

At the same time as this stitching procedure is being 
performed, catalogs containing a uniform surface den- 
sity of points placed at semi-random positions (referred 
to as 'random catalogs') are produced for the purposes 
of area calculation. These contain points placed across 
each patch at a nominal density of 0.1 arcsec"^. Where 
these fall outside the survey geometry, or within regions 
masked by artefacts, the points are removed. Thus, these 
essentially represent a sparse 'good pixel mask' for the 
survey. Future analyses requiring knowledge of the area 
considered around a given point (such as the area sur- 
veyed within X arcmin radius of a given galaxy cluster) 
can then count the number of points in the random cat- 
alog within that region to calculate the usable area. If 
higher density sampling is required, additional realisa- 
tions at the same sampling density can be generated to 
increase the resolution as needed. 

10. SUMMARY AND FUTURE WORK 

We have presented the methodology for the reduction 
and precise calibration of the second Red-sequence Clus- 



ter Survey (RCS-2). The primary purpose of this sur- 
vey is to build the largest sample of optically-selected 
galaxy clusters out to lookback times around half the 
age of the Universe. Results from cluster-finding (which 
will be presented in a forthcoming paper) show that the 
red-sequence redshifts are better than in RCS-1, due to 
the better uniformity and accuracy of the photometry in 
RCS-2, with typical {zrcs — Zspec\ ^ 0.02. A large num- 
ber of projects aimed at characterizing clusters found 
with our technique have been ongoing for several years. 
An accurate understanding of the mass-richness relation 
will allow the cluster survey to place constraints on cos- 
mological parameters. Future papers will present not 
only the cosmology results, but results on galaxy evohi- 
tion, both within the clusters and the field. In addition to 
the large number of cluster science works possible, many 
other studies including those concerned with the prop- 
erties of the Milky Way, such as the stellar populations 
and searches for dwarf galaxies, are possible. The large 
and accurate dataset provided by RCS-2 make it an im- 
portant wide-field imaging survey of the new generation. 
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